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Supernova neutrinos can easily be detected by a spherical gaseous TPC detector measuring very 
low energy nuclear recoils. The expected rates are quite large for a neutron rich target since the 
neutrino nucleus neutral current interaction yields a coherent contribution of all neutrons. As a 
matter of fact for a typical supernova at 10 kpc, about 1000 events are expected using a spherical 
detector of radius 4 m with Xe gas at a pressure of 10 Atm. A world wide network of several such 
simple, stable and low cost supernova detectors with a running time of a few centuries is quite 
feasible. 

PACS numbers: 13.15.+g, 14.60Lm, 14.60Bq, 23.40.-s, 95.55.Vj, 12.15.-y. 



I. INTRODUCTION 

Neutrinos appear to be excellent probes for studying 
the deep sky. They travel large distances with the speed 
of Hght. They can pass through obstacles, without get- 
ting distorted on their way and they are not affected by 
the presence of magnetic fields. Thus with neutrinos one 
can see much further than with light With light one can- 
not observe further than 50 Mpc (1 Mpc=3. 3x106 light 
years) . Furthermore the detection of neutrinos reveals in- 
formation about the source and more specifically about 
the source interior. Without neutrinos we would proba- 
bly know nothing about the sun's interior. Thus neutri- 
nos offer a good hope for understanding the early stages 
of a supernova. In a typical supernova an energy of about 
j^q53 gj-gg jg released in the form of neutrinos iL] , . These 
neutrinos are emitted within an interval of about 10 s af- 
ter the explosion and they travel to Earth undistorted, 
except that, on their way to Earth, they may oscillate 
into other flavors. The phenomenon of neutrino oscilla- 
tions is by now established by the observation of atmo- 
spheric neutrino oscillations Q interpreted as 
oscillations, as well as disappearance in solar neutrinos 
P|. These results have been recently confirmed by the 
KamLAND experiment ,5|, which exhibits evidence for 
reactor antineutrino disappearance. Thus for traditional 
detectors relying on the charged current interactions the 
precise event rate may depend critically on the specific 
properties of the neutrinos. The time integrated spectra 
in the case of charged current detectors, like the SNO ex- 
periment, depend on the neutrino oscillations j^. This, 
of course, may turn into an advantage for the study of 
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the neutrino properties |3 ■ An additional problem is the 
fact that the charged current cross sections depend on 
the details of the structure of the nuclei involved. 

During the last years various detectors aiming at de- 
tecting recoiling nuclei have been developed in connection 
with dark matter searches with thresholds in the 10 
keV region. Recently, however, it has become feasible to 
detect neutrinos by measuring the recoiling nucleus and 
employing gaseous detectors with much lower threshold 
energies. Thus one is able to explore the advantages of- 
fered by the neutral current interaction, exploring ideas 
put forward more than a decade ago 9J. This way the 
deduced neutrino fluxes do not depend on the neutrino 
oscillation parameters (e.g. the mixing angles). Even 
in our case, however, the obtained rates depend on the 
assumed characteristic temperature for each flavor, see 
sec. IVII . Furthermore this interaction, through its vec- 
tor component, can lead to coherence, i.e. an additive 
contribution of all nucleons in the nucleus. Since the 
vector contribution of the protons is tiny, the coherence 
is mainly due to the neutrons of the nucleus. 

In this paper we will derive the amplitude for the differ- 
ential neutrino nucleus coherent cross section. Then we 
will utilize the available information regarding the en- 
ergy spectrum of supernova neutrinos and estimate the 
expected number of events for all the noble gas targets. 
We will show that these results can be exploited by a net- 
work of small and relatively cheap spherical TPC detec- 
tors placed in various parts of the world (for a description 
of the apparatus see our earlier work pJ]||). The opera- 
tion of such devices as a network will minimize the back- 
ground problems. There is no need to go underground, 
but one may have to go sufficiently deep underwater to 
balance the high pressure of the gas target. Other types 
of detectors have also been proposed [ill, [13 ■ 
Large gaseous volumes are easily obtained by employing 
long drift technology (i.e TPC) that can provide massive 
targets by increasing the gas pressure. Combined with 
an adequate amplifying structure and low energy thresh- 
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olds, a three-dimensional reconstruction of the recoiling 
particle, electron or nucleus, can be obtained. The use 
of new micropattern detectors and especially the novel 
Micromegas [13 provide excellent spatial and time accu- 
racy that is a precious tool for pattern recognition and 
background rejection The virtue of using such 

large gaseous volumes and the new high precision mi- 
crostrip gaseous detectors has been recently discussed in 
a dedicated workshop 0| and their relevance for low 
energy neutrino physics and dark matter detection has 
been widely recognized. Such low-background low-energy 
threshold systems are actually successfully used in the 
CAST [T^ . the solar axion experiment, and are under de- 
velopment for several low energy neutrino or dark matter 
projects [13, Hi. 



II. THE NOSTOS DETECTOR NETWORK 

Before we embark on our calculations involving the 
event rates for supernova neutrino detection of a gaseous 
TPC detector we like to spend a little time in discussing 
the detector. A description of the NOSTOS project and 
details of the spherical TPC detector are given in ^3] 
and is shown inQ] We have built a spherical prototype 




Detector + tritium source 



FIG. 1: A sketch of the NOSTOS detector drawn originally for 
detecting very low energy electron recoils. The corresponding 
one for detection of nuclear recoils is analogous. 

1.3 m in diameter which is described in (see Fig. 
121) . The outer vessel is made of pure Cu (6 mm thick) 




FIG. 2; First prototype - The SACLAY sphere: R = 1.3 m, 
V — 1 m^, spherical vessel of Cu 6 mm thick, pressure up to 
5 bar (already tested up to 1.5 bar), vacuum tight, 10~^ mbar 
(outgassing ~ 10^ mbar/s) 



allowing to sustain pressures up to 5 bar. The inner de- 
tector is just a small sphere, 10 mm in diameter, made 
of stainless steel as a proportional counter located at the 
center of curvature of the TPC. We intend to use as am- 
plifying structure a spherical TPC 113 and developments 
are currently under way to build a spherical TPC detec- 
tor using new technologies. First tests were performed 
by filling the volume with argon mixtures and are quite 
promising. High gains are easily obtained and the sig- 
nal to noise is large enough for sub-keV threshold. The 
whole system looks stable and robust. The advantages 
of using the spherical detector concept are the following; 

1. The natural radial focusing of the TPC allows to 
collect and amplify the deposited charges by a sim- 
ple and robust detector using a single electronic 
channel to read out a large gaseous volume. The 
small size associated to small detector capacitance 
permits one to achieve very low electronic noise. 
In the present prototype the noise is as low as 
a few hundred electrons and has easily been ob- 
tained; with optimized low noise amplifiers we hope 
to lower it to the level of a few tenths. This is a key 
point for the obtaining a very low energy threshold, 
i.e. down to 100 eV, by operating the detector at 
moderate gain of about 100. Such low gains are 
easily obtained at atmospheric pressure and open 
the way to operate the TPC at high pressures. We 
target to achieve pressures as high as 10 bar for 
Xenon gas. Even higher pressures by a factor 3-6 
are aimed at in the case of Argon gas in order to 
achieve, to first order, the same number of events. 

2. The radial electric field, inversely proportional to 
the square of the radius, is a crucial point for mea- 
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suring the depth of the interaction by a simple anal- 
ysis of the time structure of the detector signal. A 
position resolution of about 10 cm has been already 
obtained, a fact that is of paramount importance 
for improving the time resolution of the detector 
and rejecting background events by applying fidu- 
cial cuts. 

3. Building a high pressure metallic sphere, for in- 
stance made out of stainless steel or copper, seems 
to assure an excellent quality of the gas mixture 
and turns out that a single gas filling with pure gas 
is sufhcient to maintain the stability of the signal 
for several months. We are pushing the technology 
to improve the properties of the various elements 
in order to achieve stability over many years. 

4. Big high pressure-secure tanks are under develop- 
ment by many international companies for hydro- 
gen or oil storage, and therefore the main element 
of the TPC could be shipped at moderate cost. 

Our idea is then to build several such low cost and ro- 
bust detectors and install them in several places over the 
world. First estimations show that the required back- 
ground level is modest and therefore there is no need for 
a deep underground laboratory. A mere 100 meter wa- 
ter equivalent coverage seems to be sufficient to reduce 
the cosmic muon flux at the required level (in the case 
of many such detectors in coincidence, a modest shield is 
sufficient) . The maintenance of such system could be eas- 
ily assured by Universities or even by secondary schools. 
Thanks to the simplicity of the system it could be oper- 
ated by young students with a specific running program 
and simple maintenance every a few years. Notice that 
such detector scheme, measuring low energy nuclear re- 
coils from neutrino nucleus elastic scattering, do not de- 
termine the incident neutrino vector and, therefore, it is 
not possible this way to localize the Supernova. A cluster 
of such detectors in coincidence, however, could localize 
the star by a triangulation technique. 
A network of such detectors in coincidence with a sub- 
keV threshold could also be used o observe unexpected 
low energy events. This low energy range has never been 
explored using massive detectors. A challenge of great 
importance will be the synchronization of such a detec- 
tor cluster with the astronomical 7-ray burst telescopes 
to establish whether low energy recoils are emitted in 
coincidence with the mysterious 7 bursts. 

III. STANDARD AND NON STANDARD 
WEAK INTERACTION 

The standard neutral current left handed weak inter- 
action can be cast in the form: 

£u = [vo.r{l - [ql^.{9v{q) - 9A{q)l'')q] 

(1) 



(diagonal in flavor space) with 

9v{u) = ]^-^siT?ew , 5a(u) = ^ ; 

12 1 
9v{d) = -- + -sm^ew , 9a{u) = -- (2) 

At the nuclcon level we get: 

(3) 

with 

9v{p) = ^-2siB^ew , 5a(p) = 1-27^ ; 

9vin) = -\ , 9A{n) = -1.27^ (4) 

Beyond the standard level one has further interactions 
which need not be diagonal in flavor space. Thus 

9v{q) - 9A{q)l^ 

{9r{q) - 9r{q)i') + (a'^<5„, + ef,) (1 - 7^) 

[D^ril - r'W\ -> [Pa7'^(l - (5) 
Furthermore at the nucleon level 
9y {N) - SA (iV)7^ - (.gf ^ (N) - gf ' iN)^') 6^0 

[^o.r{i-i>o]~^ [vo.r{i-i^)vp\ (6) 

with 

In the above expressions A''^ can arise, e.g., from radia- 
tive corrections, see e.g. PDG 0| and e^^ from R-parity 
violating interactions in supersymmetric models poj-[2lj|. 
Indeed since R-parity conservation has no robust theoret- 
ical motivation one may accept an extended framework 
of the MSSM with R-parity non-conservation MSSM. In 
this case the superpotential W acquires additional R- 
paity violating terms: 

+ K^^^D'^Dl+ii.L.H^ (8) 

Of interest to us here is the Xijk^iQjDk involving 
first generation quarks and s-quarks, i.e the term 
^aii^aQiDI. From this term in four component no- 
tation we get the contribution 

Kii {d'j^VaL - UROiL) d", a = e,/i, r 

where i/aL = 5(1 ~ 75)^^0 etc. Thus 



4 



The first term at tree level yields the interaction 



X' V 



2 



(9) 



By performing a Fierz transformation we can 
rewrite it in the form: 



1 ^all^'/311 



(10) 



The previous equation can be cast in the form: 



G 



F d 



V2 



^11 T 



(11) 



There is no such term associated with the u quark, 

• Proceeding in an analogous fashion the collabora- 
tive effect of the first and second term, for a, (3 — 
e, fi, r, yields the charged current contribution: 

Cdu = ^eip[ar{l--/')M [u7M(l-7')rf](12) 

• Finally the second term, for a, (3 = e, fj,, r, leads to a 
neutral current contribution of the charged leptons: 

£„ = ^4^[aj'^il-j^)(3][u^,{l-j')u] (13) 

The above non standard flavor changing neutral current 
interaction have been found to play an important role in 
the in the infall stage of a stellar collapse ,2^]. Further- 
more precise measurements involving the neutral current 
neutrino-nucleus interactions may yield valuable infor- 
mation about the non standard interactions 24] . They 
are not, however, going to be further considered in this 
work. 



momentum transfer to the nucleon they take the form 
(see previous section): 



5y = -2 sin^ 6*,^ + 1/2 w 0.04 , gA 



1.27 



iiy,p) (15) 



gv = -1/2 , gA 



1.27 



(z.,n) (16) 



In the above expressions for the axial current the renor- 
malization in going from the quark to the nucleon level 
was taken into account. For antincutrinos gA ~9A- 
To set the scale we write: 



^ 5.14 X 10-^1 ^ 
27r MeV 



(17) 



The nucleon energy depends on the neutrino energy and 
the scattering angle, the angle between the direction of 
the recoiling particle and that of the incident neutrino. 
In the laboratory frame it is given by: 



Tn = 



2 'mjM{Ey cos 6*)^ 



[niN + E^Y-{E^ cos 



, < 61 < 7r/2 (18) 



(forward scattering). For sufficiently small neutrino en- 
ergies, the last equation can be simplified as follows: 



Tn 



2{E^cos0f 
ruN 



The above formula can be generalized to any target and 
can be written in dimensionless form as follows: 



2cos^ 



y = 



(1 + l/x^)^ - C0S2 I 



, y 



recoil 



E, 



^recoil ^recoil 

(19) 

The maximum energy occurs when 9 = 0, Umax = 
(i+i/x y^~i ' ™ agreement with Eq. (2.5) of ref. [1|. This 
relationship is plotted in Fig. El One can invert Eq. El 



IV. ELASTIC NEUTRINO NUCLEON 
SCATTERING 

The cross section for elastic neutrino nucleon scattering 
has extensively been studied. It has been shown that at 
low energies the weak differential cross section can be 
simplified and be cast in the form: (li|.|i22l]: 



da 
dTN 



2tt 

i9v-9Ani 



[{gv + 9A? 



(14) 



N^2 



E, 



{g\ 



2 \ ^nTn 1 



g'v) 



El 



where rriAr is the nucleon mass and gv, gA are the weak 
coupling constants. Neglecting their dependence on the 
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FIG. 3: The maximum recoil energy as a function of the neu- 
trino energy (both in units of the recoiling mass). The scale 
is realistic for nuclear recoils of the type of experiments con- 
sidered here. 
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and get the neutrino energy associated with a given recoil 
energy and scattering angle. One finds 



-1 + cos 61 



2 

1 + - 

V\ 



< 61 < 7r/2 (20) 



The minimum neutrino energy for a give recoil energy is 
given by: 



-1 



1 + ^ 
V 



2^ 



1 + ^) 

V 



(21) 



in agreement with Eq. (4.2) of ref. The last equa- 
tion is useful in obtaining the differential cross section 
(with respect to the recoil energy) after folding with the 
neutrino spectrum and it is shown in Fig. 0] 




FIG. 4: The minimum neutrino energy as a function of tlie 
recoil energy (botli in units of the recoiling mass). Tlie scale is 
not realistic for nuclear recoils considered here, but relevant 
for other experiments. The realistic scale for experiments 
considered here can be deduced from Fig. 0|by exchanging 
the coordinate axes. 



V. COHERENT NEUTRINO NUCLEUS 
SCATTERING 

From the above expressions we see that the vector cur- 
rent contribution, which may lead to coherence, is negli- 
gible in the case of the protons. Thus the coherent contri- 
bution [25 may come from the neutrons and is expected 
to be proportional to the square of the neutron number. 
The neutrino- nucleus scattering can be obtained from the 
amplitude of the neutrino nucleon scattering under the 
following assumptions: 



• Employ the appropriate kinematics, i.e. 
volving the elastically scattered nucleus. 



those in- 



• Ignore effects of the nuclear form factor. Such ef- 
fects, which are not expected to be very large, are 
currently under study and they will appear else- 
where. 



• The effective neutrino-nucleon amplitude is ob- 
tained as above with the substitution 



P 

A 



E 



N 



"N 



Ea 
A 



with q the nucleon momentum and p the nuclear 
momentum. 

Under the above assumptions the neutrino-nucleus cross 
section takes the form: 



+ {Mv-MAf{l 



Gl^AniN 
X 

[(Mv + MaY 



1 



A-ITa 



A E, 
A-ITa 1 



+ {Mi- Mi.) 



A 771 jv E^/Ta 
2 - AtumTa , 



El 



- 1 
(22) 



Where My and Ma are the nuclear matrix elements asso- 
ciated with the vector and the axial current respectively 
and Ta is the energy of the recoiling nucleus. The axial 
current contribution vanishes for 0^ =J> 0+ transitions. 
Anyway it is negligible in front of the coherent scattering 
due to neutrons. Thus the previous formula is reduced 
to: 



\dTAj, 



G%AmN 



2tt 



iNyA)F,ohiA,TA,E,), 



F,,hiA,TA,E,) 
A-ITa 



A-ITa 



I 



1 



A E^J ' E^ 
AtiinTa 



A ruN E^/Ta - 1 



El 



(23) 



The function Fcoh{A, Ta, E^) is shown in FigElas a func- 
tion of the recoil energy in the case of Ar (N=22) and Xe 
(N=77) for 10, 20, 30 and 40 MeV respectively. We see 
two reasons for enhancement of the cross section: 

• The overall A factor due to the kinematics, which 
is counteracted by the smaller nuclear recoil en- 
ergy when compared to the nucleon recoil energy 
for the same neutrino energy. This factor will be 
absorbed into the energy integrals, see the function 
Ffoid{A,T, {TaU) below. 

• The iV^ enhancement due to coherence. 



VI. 



SUPERNOVA NEUTRINOS 



The number of neutrino events for a given detector 
depends on the neutrino spectrum and the distance of 
the source. We will consider a typical case of a source 
which is about 10 kpc. I.e. D = 3.1 x 10^^ cm ( of the 
order of the radius of the galaxy) with an energy output 
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FIG. 5: The function Fcoh{A,TA, E^) as a function of the 
recoil energy Ta for, from left to right, Ei, = 10, 20, 30, 40 
MeV. The results shown are for Xe on the top and Ar at the 
bottom 



of 3 X 10^^ ergs with a duration of about 10 s. 
The neutrino spectra are parametrized as follows: 



dN 



/3f 



E, 



Exp -A 



E, 



E^ y 
(24) 

where -< Ei >- is the average energy of neutrino flavor i 
with E^^ < Ep^ < E,y^{Ep^) for x — fi,T. The param- 
eters $i, -< Ei >- and Pi for each flavor, Ve^^e, Vx (j^x ) 
with X = fi,T are determined phenomenologically '2^. 
In the present paper, in order to minimize the number 
of parameters, we will assume for simplicity that each 
neutrino flavor is characterized by a Fermi-Dirac like dis- 
tribution times its characteristic cross section, which is 
adequate for our purposes. Thus: 



^2 



A 



dN _ 

dK~'^ 1 + exp{E^/T) ^ JTl + e^ 



= 

^ T 
(25) 

with J = ^§2^ ^ ^ constant and T the temperature of 
the emitted neutrino flavor. Each flavor is characterized 
by its own temperature as follows: 

T = 8 MeV for v^,, Vr,Vt,, Vr,T = 5 (3.5) MeV for Ve {v^) 

The constant A is determined by the requirement that 
the distribution yields the total energy of each neutrino 
species. 



Uy = — I dx ^ 



1 -I- T 



We will further assume that C/^ = 0.5 x 10^"^ ergs per 
neutrino flavor. Thus one flnds: 

A = 0.89x10^^** (i^e), 0.63x10^^ (z>e) ,0.39x10^^ (other) 
The emitted neutrino spectrum is shown in Fig. The 




20 40 

E^ MeV 

FIG. 6: The supernova neutrino spectrum. The short dash, 
long dash and continuous curve correspond to i/ejt'c and all 
other flavors respectively 



differential event rate (with respect to the recoil energy) 
is proportional to the quantity: 



dR _ \{T) 
dJ^ ~ J 



/ dxFcohiATA,xT) — 
Jo ^ 



(26) 



with A(r) = (0.89,0.63,0.39) for i^e,i^e and aU other fla- 
vors respectively. This is shown in Figs. |7|and[71 The 
total number of expected events for each neutrino species 
can be cast in the form: 



No of events 



HAT, (T. 



A)th) 



C,{T)h{A, T, (TaU), 

Ffoid{A,T, {TaU) 
FfoM{^Q,T, (TaU) 



(27) 
(28) 



with 



FfoidiAT, [Ta. 



tH) - -J 



{Ta 



dTA 
IMeV 



dxFcoh{A,TA,xT) 



1 



■(29) 



and 



CAT) 



G^ruNlMeV 

'hi r 



A(r) 



1 PV 



(30) 



Where k is Boltzmann's constant, P the pressure, V the 
volume, and Tq the temperature of the gas. 

Summing over all the neutrino species we can write: 

Noofevents = ar(A)f^i^^Qu(A) (31) 
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with 



— > McV 



FIG. 7: The differential event rate as a function of the recoil 
energy Ta, in arbitrary units, for Xe. On the top we show the 
results without quenching, while at the bottom the quenching 
factor is included. We notice that the effect of quenching is 
more prevalent at low energies. The notation for each neu- 
trino species is the same as in Fig. |S| 




McV 
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FIG. 8: The same as in Fig. |7|for the Ar target. 



' N 



0.5 X lO^^ergs \ D ) 10 Atm 



R 

4m 



300 



(32) 



In the above expression r{A) is a kinematical parameter 
depending on the nuclear mass number, which is essen- 
tially unity. 

K{A, {TA)th) is the rate at a given threshold energy 
divided by that at zero threshold. It depends on the 
threshold energy, the assumed quenching factor and the 
nuclear mass number. It is unity at (T^)t/i) = 0. The 
function r{A) is plotted in |51 It is seen that it can be 
well approximated by unity. 

From the above equation we find that, ignoring quench- 
ing, the following expected number of events: 

1.25, 31.6, 153, 614, 1880 for He, Ne, Ar, Kr and Xe 

(33) 

respectively. For other possible targets the rates can be 
found by the above formulas or interpolation. 
The function K{A, {TA)th) is plotted in Fig. UHl for 
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FIG. 9: The function r{A) versus the nuclear mass number. 
To a good approximation r{A) ~ 1.0 (for definitions see text) 



threshold energies up to 2keV. We see that the threshold 
effects are stronger in heavier systems since, on the aver- 
age, the transfered energy is smaller. Thus for a threshold 
energy of 2 keV in the case of Xe the number of events 
is reduced by 30% compared to those at zero threshold. 
The quantity Qu{A) is a factor less than one multiply- 
ing the total rate, assuming a threshold energy {TA)th — 
lOOcV, due to the quenching. The idea of quenching is 
introduced, since, for low emery recoils, only a fraction 
of the total deposited energy goes into ionization. The 
ratio of the amount of ionization induced in the gas due 
to nuclear recoil to the amount of ionization induced by 
an electron of the same kinetic energy is referred to as 
a quenching factor Qfac- This factor depends mainly on 
the detector material, the recoiling energy as well as the 
process considered [2g. In our estimate of Qu{Ta) we as- 
sumed a quenching factor of the following empirical form 
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0.0005 0.001 0.0015 0.002 

[TaU ^MeV 

FIG. 10: The function K{A, {TA)th) versus {TA)th for various 
nuclear mass numbers without the quenching factor. From 
top to bottom He, Ne, Ar, Kr and Xe. (for definitions see 
text) 
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motivated by the Lidhard theory |2(i|-|28|: 

" Ta 



Qfac{TA) = n 



IkeV 



n ~ 0.256 , r2 ~ 0.153 

(34) 



Then the parameter Qu{A) takes the values: 

0.49, 0.38, 0.35, 0.31, 0.29 for He, Ne, Ar, Kr and Xe 

(35) 

respectively. The effect of quenching is larger in the case 
of heavy targets, since, for a given neutrino energy, the 
energy of the recoiling nucleus is smaller. Thus the num- 
ber of expected events for Xe assuming a threshold energy 
of 100 eV is reduced to about 560. 

The effect of quenching is exhibited in Fig ^] for the two 
interesting targets Ar and Xe. 

We should mention that it is of paramount importance 
to experimentally measure the quenching factor. The 
above estimates were based on the assumption of a pure 
gas. In our detection scheme the Xe gas carrier (A) is 
mixed with a small fraction of low ionization potential 
gas (B). Thus a part of the excitations produced on the 
Xe atoms could be transferred to ionization through the 
well known Penning effect as follows: 



A* +B 



A + B* 



(36) 



Such an effect will lead to an increase in the quenching 
factor and needs be measured. 



VII. CONCLUSIONS 

In the present study it has been shown that it is 
quite simple to detect typical supernova neutrinos in our 
galaxy, provided that such a supernova explosion takes 
place (one explosion every 30 years is estimated [31^). 
The idea is to employ a small size spherical TPC detector 
filled with a high pressure noble gas. An enhancement of 



b 0.75 - 
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FIG. 11: The function K{A, {TA)th) versus (Ta)*^ for the 
target Ar on the top and Xe at the bottom. The short and 
long dash correspond to no quenching and quenching factor 
respectively. One sees that the effect of quenching is less pro- 
nounced at higher thresholds. The differences appear small, 
since we present here only the ratio of the rates to that at zero 
threshold. The effect of quenching at some specific threshold 
energy is not shown here. For a threshold energy of 100 eV 
the rates are quenched by factors of 3 and 3.5 for Ar and Xe 
respectively (see Eq. I35II . 



the neutral current component is achieved via the coher- 
ent effect of all neutrons in the target. Thus employing, 
e.g., Xe at 10 Atm, with a feasible threshold energy of 
about 100 eV in the detection the recoiling nuclei, one 
expects between 600 and 1900 events, depending on the 
quenching factor. We believe that networks of such ded- 
icated detectors, made out of simple, robust and cheap 
technology, can be simply managed by an international 
scientific consortium and operated by students. This net- 
work comprises a system, which can be maintained for 
several decades (or even centuries) . This is is a key point 
towards being able to observe few galactic supernova ex- 
plosions. 
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